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Abstract 

In this paper we study the Wiener-Hopf factorization for a class of Levy processes with double- 
sided jumps, characterized by the fact that the density of the Levy measure is given by an infinite 
series of exponential functions with positive coefficients. We express the Wiener-Hopf factors as 
infinite products over roots of a certain transcendental equation, and provide a series representation 
for the distribution of the supremum/infimum process evaluated at an independent exponential 
time. We also introduce five eight-parameter families of Levy processes, defined by the fact that 
the density of the Levy measure is a (fractional) derivative of the theta-function, and we show 
that these processes can have a wide range of behavior of small jumps. These families of processes 
are of particular interest for applications, since the characteristic exponent has a simple expression, 
which allows efficient numerical computation of the Wiener-Hopf factors and distributions of various 
functionals of the process. 
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1 Introduction 



Wiener-Hopf factorizaton and related fluctuation idenitities (see [4], [9], [12] and [18]) allow us to study 
various functionals of a Levy process, such as extrema, first passage time, overshoots and undershoots, 
etc. There is a growing number of applications of these functionals and of the Wiener-Hopf factorization 
techniques in many areas of Applied Probability, most prominently in Mathematical Finance and Insur- 
ance Mathematics. However the number of processes with two-sided jumps for which one can obtain 
explicit results on fluctuation identities is by no means large. In fact until very recent times the only 
known examples consisted of a dense subclass of stable processes (see [8]), processes with phase-type 
distributed jumps (see [1] or [2] ), and more generally, processes having positive jumps with rational 
transform and arbitrary negative jumps (see [16]). In the last several years there have appeared a num- 
ber of new results on the Wiener-Hopf factorization with varying degree of explicitness. First of all 
we would like to mention the class of Lamperti-stable processes, which can be obtained by the Lam- 
perti transformation (see [14]) from positive self-similar Markov processes, which are related to stable 
processes. The family of Lamperti-stable processes is very interesting (and quite unique), since many 
fluctuation identities can be obtained in closed form in terms of elementary or special functions, see [5], 
[6], [13] and the references therein. Second, in [11] we have introduced a ten-parameter family of Levy 
processes, for which the Wiener-Hopf factors can be identified in a semi-explicit form, more precisely as 
infinite products involving solutions to a certain transcendental equation. It is the goal of this paper 
to build on our earlier ideas and to extend the results in [11], as well as to provide new analytically 
tractable examples of Levy processes with semi-explicit Wiener-Hopf factorization. 

The main idea in [11] consists in the following observation: if the characteristic exponent "^{z) can be 
extended to a meromorphic function (which means that its only singularities in the complex plane are 
poles) and if we have some additional information about the asymptotic behavior of the solutions of the 
equation ^(^) + q = 0, then we can obtain the Wiener-Hopf factors in essentially the same way as if "^{z) 
was a rational function. Our first contribution in this paper is to show that the same technique works 
in a much more general setting: if the Levy measure is an infinite series of exponential functions with 
positive coefficients then the Wiener-Hopf factorization and the distribution of extrema can be computed 
rather explicitly in terms of the solutions of the equation \l/(z) + g = 0. As we will see, the following 
property plays the most important role in the proof of our main result: the Levy measure is an infinite 
series of exponential functions with positive coefficients if and only if for every g > the zeros and poles 
of "^{z) + q lie on the imaginary axis and interlace. As our second contribution we present several new 
families of Levy processes which satisfy the above condition and which will be interesting for applications 
due to their analytic tractability. 

The paper is organized as follows: in section 2 we prove Theorem 1 on the Wiener-Hopf factoriza- 
tion and distribution of extrema, while in section 3 we construct five eight-parameter families of Levy 
processes, for which the density of the Levy measure is essentially the fractional derivative of the theta 
function. For each of these five families we provide explicit formulas for the characteristic exponent 
and asymptotic expressions for the large roots of g + "^{z) = 0, which are important for the efficient 
implementation of numerical algorithms. 
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2 Main Results 

First let us present several definitions and notations which will be used in this paper. We define the 
open/closed upper half plane as 

C+ = {2 e C : lm{z) > 0}, C+ = {2 G C : lm{z) > 0}, 

and similarly for the negative half plane and for the positive/negative real half line. We will study a one- 
dimensional Levy process X started from zero, which is defined by the characteristic triple {fi, a, 11). The 
characteristic exponent "^{z) = — ln(E[exp(i2;Xi)]) can be computed via the Levy-Khintchine formula 
(see [4]) as folios 

^(2) = ^ah^ -ifiz- j (e'-'^' - 1 - izxh{x)) U{dx), zeR, (1) 

where h{x) is the cutoff function (everywhere in this paper we will use h{x) = or h{x) = 1 as the 
measure n(da;) will have exponential tails). We define the supremum/infimum processes as 

St = sup Xs, It = inf Xs. 

0<s<t 0<s<t 

Wiener- Hopf factors (see [4], [9], [12] or [18]) are defined as 

(f)+{z) = E [e^^^-w] , for 2 e C+, 0^ (^) = E [e^^^-w] , for ;2 e C", 

where the random variable e(g) is exponentially distributed with parameter g > and is independent of 
the process X. 

In order to specify the Levy measure n(da;) we start with the four sequences of positive numbers 
{an, Pn, o,n, Pn}n>i, and assume that the sequences {pn}n>i and {p„}„>i are stricly increasing, and that 
Pn — )■ +00 and Pn — +00 as n ^ +00. 

Assumption 1. Series J2n>i^nPn^ f*^'^ Xln>i ^"Pn^ converge. 
Next we define the function 7t{x) as 

7r(x) = I{x >0)J2 (^nPne-'-'^ + l{x <0)Y, «nPne^"^ (2) 

n>l n>l 

It is easy to see that the above series converges for all x ^ 0, uniformly in a; on M \ (— e, e) for all e > 0. 
To check this one should use Assumption 1 and the inequality exp(— p„e) < which holds for all n 
large enough. It is clear that n{x) is a positive function which decays exponentially as |x| — )■ oo, and 
as we show in the next Proposition, Assumption 1 guarantees that 7r(x) can be used to define a Levy 
measure. 

Proposition 1. Assumption 1 implies that j^x'^'K{x)(\x < oo. 

Proof. As we have mentioned above, for all e > the first series in (2) converges uniformly for x G (e, oo), 
thus we can integrate term by term and obtain 



oo 



/x^7r(x)dx = ^^flnPn / ''"^^x = a„p„ ^ / u^e "dw 
'J ^"^1 'J 



n>l n>\ 
— e - PnC 



oo 

< XI ' / ^'e-"du < 2 X anPl 

n>\ y n>l 
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where in the second step we have changed the variable of integration x i— >■ m = p„x. We see that the 
integral on the left-hand side of the above inequality increases and is bounded as e — )■ O"*", therefore it 
converges. The convergence of the integral over (— cxo, 0) is proved similarly. □ 

Proposition 1 allows us to use tt{x) as the density of a Levy measure n(dx) = 7r(x)da;. Note that 
since 7r(x) decays exponentially as x — t- oo we can use the cutoff function h{x) = 1 in (1); this will be the 
default choice everywhere in this section. The above property also allows us to work with the Laplace 
exponent, defined as 

= In (E [e'^']) = -^{-iz). (3) 

Proposition 2. The Laplace exponent (j){z) is a real meromorphic function which has the following 
partial fraction decomposition 

= V + + ^ , +z'Y. ^r\v (4) 

2 Pn{Pn -Z) ^ Pn[Pn + z) 

Proof. Use (1), (2) and (3), integrate term by term and rearrange the resulting infinite series. □ 

Proposition 3. Assume that g > 0. Equation (j)(z) = q has solutions {(n, —Cn}n>i, where {Cn}n>i O'lT'd 
{Cn}n>i sequences of positive numbers which satisfy the following interlacing property 

< Ci < Pi < C2 < P2 < • • • (5) 
< Ci < pi < C2 < P2 < • • • 

Proof. Using (4) we rewrite the equation (f){z) = g as 

Z'Y Qn 1^2^2_ /g) 

^^PniPn-z) ^^PniPn + z) ^ 2 ^ 

Denote the left-hand side of (6) as (piiz) and the right-hand side as (f)2{z)- First let us check that there 
exists a solution to (6) on intervals (0, pi) and (— pi, 0). We observe that 01 (0) = and (pi{z) /• +oo as 
z pi or as 2; \ — pi; at the same time function </'2(;z) is continuous and 02(0) > 0, all that is left to do 
is to apply the Intermediate Value Theorem. Other intervals can be verified in a similar way. □ 

It is important to emphasize that Proposition 3 does not claim that the sequences {Cn, — Cn}n>i 
include all solutions to 0(z) = q. The statement is that some solutions are real and that they interlace 
with the poles {p„, — p„,}n>i- As we will see later, it is in fact true that there are no other solutions, but 
the proof of this statement is not trivial and requires some deep results from the theory of meromorphic 
functions. . 

Our first main result in this paper is the following Theorem, which identifies the Wiener-Hopf factors 
and the distribution of extrema of the process X. 

Theorem 1. Assume that g > 0. Then for Re(z) > 

n>l C" n>l C,n 
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The distribution of S^^g) can be identified as an infinite mixture of exponential distributions 



P(5e(,) = 0) = Co, ^P('5e(g) <X) = J2 CnCne^^"", X > 0, 



(8) 



n>l 



where the coefficients {c„}„>o are positive, satisfy J2n>o^"- ~ ^' ^'^^ computed as 



Co = lim TT — , Cn 

n^^r^ J. J. 



k=l 



Cn 
Pn 



k>l Cfe 



(9) 



r/ie distribution of —Ic(q) has the same form as above, with {pnXn} replaced by {pn,Cn}- 

Our proof of the Theorem 1 is based on the two important results from the theory of meromorphic 
functions, which we present below. These two theorems combined together state that the Levy measure 
has the form (2) if and only if for every q > the zeros and poles of the real meromorphic function 
~ Q rsal and interlace, and this function can be represented as an infinite product. 

Theorem 2. ([15], page 220, Theorem 1) A real meromorphic function f{z) satisfies f{z) G C+ for all 
z E if and only if it can be represented in the form 



f{z) = c 



z — ao 
z -bo 



1 



n 



(10) 



^6Z\{0} 

where c > 0, the zeros a„ and the poles bn are real, satisfy the interlacing property 

bn < an < bn+i, n G Z, 

and a_i < < 6i. 

Theorem 3. ([7], page 197, Theorem 1) A real meromorphic function f{z) satisfies f{z) G C+ for all 
z E'C'^ if and only if it can be represented in the form 



f{z) = az + P + -^+ V Bn 
bn - z ^ 



nGZ\{0} 



1 



6„ - z b„ 



where the poles bn are real, bn < bn+i and 6_i < < 6i, a > 0, /3 G M, Bn > for n > and the series 
^ Bnbn^ converges. 

Proof of Theorem 1: First let us prove that the infinite products in (7) converge. The product Yl^n 
converges if and only if the series Yli^n — 1) converges, thus we need to consider the series 



E 

n>l 



pn 



1 + f 



1 1 



(12) 



Using the interlacing property (5) we find that 

< ■ ■ ■ < P2 ^ < C2^^ < pr^ < cr\ 
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thus we obtain 



^ VCn PnJ ^ VCn Cn+1 / 
n>l ^ ' n>l ^ 3 1/ 



and conclude that the series in the right-hand side of (12) converges, which guarantees the convergence 
of the infinite products in (7). 

Next, let us establish the following infinite product factorization 

q „ i - - „ i 



^ ' n>l Cn n>l C„ 

To prove this factorization, we use (4) and rewrite the function (0(2;) — in the following form 
= -(T2; + /i h> a„ + > a„ — ^ ^ . 

Using the above equation and Theorem 3 we conclude that function (0(2) — q)jz maps the upper half 
plane into itself, therefore applying Theorem 2 and Proposition 3 we find that 



{z)-q z-C 



1 - ^ „ 1 



n>l Pn n>l Pn 



which in turn implies (13). 
Let us introduce functions 



/(^) =111^4' ^»=nT^- (14) 

n>l Cii n>l 

Using Theorem 2 we find that the real meromorphic function zf{z) maps the upper halfplane into itself, 
thus again we apply Theorem 3 to obtain that zf{z) has the partial fraction decomposition (11), which 
implies that f{z) has the partial fraction decomposition 

n>l '"^ + ^ 

where c„ > for > and the series ^ c„ converges. Setting 2; = in the above equation we find 
that ^ c„ = 1, and therefore f{z) is the Laplace transform of the mixture of exponential distributions, 
which implies that f{z) is the Laplace transform of a positive infinitely divisible random variable with 
zero drift. The same result holds for f{z). Equation (13) tells us that q/{q — 4>{z)) = f{—z)f{z), and 
using the uniqueness of the Wiener- Hopf factorization (see [18]) we conclude that (f)g{iz) = f{z) and 
(f)g{—iz) = f{z), thus obtaining (7). 

In order to finish the proof we only need to establish formulas (9). The expression for cq follows from 
the following computations 

^ 1 _|_ ^ TV . 

= 0) = lim E \e-'^<^^] = lim lim TT ^ = lim TT 

L 7V-i.+oo J- J- 1 + ^ 7V-i.+oo J- J- 

n=l 



',n n=l 
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where in the last step we have interchanged the two hmits, which is possible due to the fact that the 
limit 

lim TT^^ 

71=1 t,n 

converges uniformly in N. The expression for c„ in (9) follows from (14) and the fact that c„^„ is the 
residue of f{z) at the pole z = —(n- □ 



3 Examples 

Theorem 1 gives us the Wiener-Hopf factors and the distribution of Se{q) and Ie{q) for any Levy process 
whose Levy measure is an infinite series of exponential functions with positive coefficients. All the 
formulas in Theorem 1 are based only on the zeros and the poles of the meromorphic function (f){z) — q. 
The poles — p„}„>i are usually known explicitly, but in order to find the zeros {(n, —Cn}n>i one has to 
solve the transcendental equation (f){z) = q, and this has to be done numerically. It is clear that if we have 
to rely on the partial fraction decomposition (4) in order to evaluate 4>{z), such numerical computations 
will be quite challenging, if not impossible. Thus it is very important to find families of Levy processes, 
for which the Levy measure is an infinite series of exponentials, and at the same time the Laplace 
exponent can be computed in closed form. In [11] we have introduced the /3-family: a ten-parameter 
family of Levy processes, for which the characteristic exponent can be computed in terms of the beta 
function. A very special subclass of the /3-family is a five-parameter family of Levy processes (see Section 
3 in [11]), whose jump component behaves similarly to the normal inverse Gaussian process (see [3]), and 
the characteristic exponent is expressed in terms of elementary (trigonometric) functions. In this section 
we present five parametric families of Levy processes, which have a number of desirable properties: the 
characteristic exponent is expressed in terms of rather simple functions, such as trigonometric functions 
or the digamma function ip{z) = T'{z)/T{z) (see [10] for the definition and properties of the digamma 
function). The density of the Levy measure decays exponentially at infinity and has a singularity at zero 

of the order x ^ {|) 1; f; 2, |}, thus "covering" the complete range of admissible singularities x ^ (0, 3). 
In particular, when x = 1 we have a process of infinite activity of jumps but of finite variation, whose 
jump part is similar to the Variance Gamma process (see [17]), and when x = 2 we obtain a process 
with infinite variation of jumps, whose jump part is similar to the normal inverse Gaussian process. 
To define these processes we introduce the function 6fc(x) as follows 

efc(x)=4,o + 2 5^n2V"'^ x>0, (16) 

n>l 

where 5^,0 = 1 if k = and 6kfl = otherwise. Note that Qo{x) = 63(0, e~^) (see [10] for the definition 
and properties of the theta functions 9i{z; r)), and for /c G N the function Ofc(x) is just the fc-th order 
derivative of the theta function ^3(0,6"^). In fact, results similar to the ones presented in this section 
can be established if we use 6*2(0,6"^) instead of 6*3(0,6"^). 
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Definition 1. For < x < 3 and x we define 

ttJx) = I(x > 0)ci/3ie-°i"efc(x/3i) < 0)c2/32e°2"efc(-x/32), 

where Ci, ai, A > and /c = x — 1/2- 

Proposition 4. Function vr^(x) /ias the following asymptotics 

%W ~ r(x)c2/32'"''|a;|-^, X ^ 0-. 

Proof. Let /i = ^/x. Then using the fact that a definite integral is the limit of the Riemann sum we 
obtain that as /i — O"*" 



E 

n>l 



n>l 





oo 








/ 




.0 



y^^e-y'dy + o(l) 



To finish the proof use Definition 1 and the above asymptotic relation. 



□ 



Proposition 4 and Definition 1 guarantee that Jj^ x^T\'^{x)dx exists, thus n^x) can be used to define the 
density of a Levy measure n^(dx) = 7r^(x)dx. We define the Levy process X using the the characteristic 
triplet (yU, 0", H-^) and the Levy-Khinchine formula (4). 

Remark 1. Note that in a similar way one can construct Levy processes with asymmetric behavior of 
small positive/negative jumps. One should extend Definition 1 to allow parameter ki {/C2} to control the 
order of the singularity of 7r{x) as x — j- 0"^ {x — )■ 0~}. 

In the next five sections we will present the results for the Levy process X with x ^ {|; 1? |;2, |}. 
The explicit formulas for the characteristic exponent "^{z) are derived using the following approach: first 
we obtain the characteristic exponent "^{z) for X = | (this is just a simple application of formulas 6.162 
in [10]) and for x = 1, in which case we have to use the following series expansion for the digamma 
function (see formula 8.362.1 in [10]) 



ip{x) 



1 

-7 - - 

X 



E 

n>l 



n{n + x) 



All the other cases, when x ^ {§5 2, |}, can be easily derived from the above two by using Definition 
1, the fact that 0^+1 (x) = — ^0fc(x) (which follows from (16)) and applying the following Proposition, 
(which can be easily established by integration by parts). 

Proposition 5. 

(i) Assume that J^_^ e~"^7r(x)dx < 00. Then for 2; G M 



(e'"^ - 1) e-''^7i'{x)dx = {a- iz)fi{z) - a/i(0), 

R+ 

where fi{z) = e'^^~"^7r(x)da;. 
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(ii) Assume that J^^ xe "^7r(x)dx < oo. Then for 2; G M 

(e'"^ - 1 - izx) e-"^7r'(x)dx = {a - iz)f2{z) - azf^{0), 



where f2{z) = (e'^^ ~ ^) ^ "^7r(x)dx. 

Asymptotic expansions for the large solutions of the equation (f){() = q are obtained using exactly the 
same technique as in the proof of the Theorem 5 in [11], and since the derivation of these expressions is 
rather lengthy and tedious we have decided to omit it. The interested reader who decides to verify these 
formulas might want to use a symbolic computation software to make the algebraic manipulations more 
enjoyable. 

All the formulas below involve the parameter a > 0, positive numbers Pi, Ci which define the Levy 
measure 11;^ via Definition 1, and two additional parameters 7 and p. The parameter 7 is uniquely 
determined by the condition ^'(O) = and the parameter p is responsible for the linear drift. In all 
the cases when x < 2, the jump part of the process has bounded variation, thus we can take the cutoff 
function in the Levy-Khintchine formula (1) as h{x) = 0, and then we have p = p. When x > 2 we 
take h{x) = 1 and then p can be uniquely expressed in terms of the characteristic triplet (/i, a, U^) via 
condition E[Xi] = i^'(O) = /i. 

3.1 The family of processes with x = 1/2 

(i) The characteristic exponent can be computed as follows 

"^{z) = ^cr^z^ - ipz - CiTT ((ai - iz)f3{^y^ coth ^tt^/ («i - iz)/]^^^ 
- C27r ((a2 + iz)f32^) ^ coth ^vry^ (02 + iz)(32^^ 



7- 



(ii) For n > 1 we have pn = ai + Pi{n — 1)^ and Pn = ^2 + P2{n — 1)^. 

(iii) If cr 7^ then the large positive solutions to 0(C) = q satisfy 

4 C2 8 C2 

C ~ P2n^ + 0:2 + -^7rn~'^ + —r^{p - a2a'^)rr^ + O in~^) , n -> +00. 
P2 P2 

(iv) If cr = and p then the large positive solutions to (f){() = q satisfy 

2c2 2 C2 

C ~ I32n^ + 0:2 H olT^^^"^ + 7 + q)n~^ + O [n~^) , n +00. 

p p P2 

3.2 The family of processes with x = 1 

(i) The characteristic exponent can be computed as follows 

^(z) = ^cr^z^ - ipz + Clip (^i\J (tti - i^)/3r^^ + Ciip («i - i^)/?!"^ 

7- 



+ C21P (^^J (02 + iz)[32 + C2i! |^-iy^(a2 + i2)/32 - 
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(ii) For n > 1 we have p„ = ai + Piin? and /)„ = 02 + /32'^^- 

(iii) If a 7^ then the large positive solutions to 0(C) = q satisfy 

4 C2 8 C2 

C ~ I32n^ + 02 H — ^^'^"^ H — 7^(/^ - a20-^)n"^ + O (n~'^) , n +00. 
0" P2 cr P2 

(iv) If cr = and /i 7^ then the large positive solutions to </){() = q satisfy 

C ~ P2n'^ + ^2 - — n"^ + ^-^(2(ci + C2) ln(n) + Co)n"^ + O (n"^ ln(n)) , n +00, 



where Cq = /iQ!2 - 7 + g + Ci In j . 



3.3 The family of processes with x = 3/2 

(i) The characteristic exponent can be computed as follows 

= -(7^2;^ — ipz + civr-y/ (ai — iz)(3i^ coth ^^r-y/ (ai — iz)/3f 
+ 0271^(02"+^^)^^ coth |^7ry'^(a7-M^)!5^^ - 7- 

(ii) For n > 1 we have p„ = cti + /3in^ and Pn = «2 + ■ 

(iii) If o" 7^ then the large positive solutions to (j){Q = q satisfy 

4 C2 8 C2 

C ~ /32«^ + tt2 + ^^^"^ + ~ ot2<y'^)n^^ + O {n^^) , n +00. 

(T P2 ^ P2 

(iv) If cr = and /i 7^ then the large positive solutions to 0(C) = q satisfy 

2c2 27r C1C2 



3.4 The family of processes with x = 2 

(i) The characteristic exponent can be computed as follows 



^ + O , n -> +00. 



-(7^2;^ — ipz — Ci(ai — i-z)/?! ^ 



02(02 + iz)l32 ^ 



i) (isJ{a2 + iz)p:,A +i)(-isj{a2 + iz)p:,^ 



+ 7- 



(ii) For n > 1 we have p„ = ai + /3in^ and Pn = c>^2 + /32'^^- 

(iii) If (7 7^ then the large positive solutions to 0(C) = q satisfy 

C ~/32ri^ + a2 + ^T-^n^^ + (n~3ln(n)) , 



-00. 
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3.5 The family of processes with x = 5/2 



(i) The characteristic exponent can be computed as follows 

"^{z) = ^cr^^^ - ipz - ciTT ((«i - i-2)/3i"^) ^ coth (ai - iz)f3{^^ 
- C2vr ((a2 + i2;)/32'^) ^ coth f 7ry'^(a7-M^)^^ ) + 7- 



(ii) For n > 1 we have pn = ai + Pin^ and /)„ = 02 + (32^^- 

(iii) If cr 7^ then the large positive solutions to 0(C) = q satisfy 

. „ 2 4 C2 Stt C1C2 _i „ / _2\ 

C~/32^ +«2 + -^-5---i 3 ^ +0{n ), n^+oo. 

P2 cr (/3i/32)2 

(iv) If cr = then the large positive solutions to (f){Q = q satisfy 

C~/32(n + ^/;o)' + «2 + ^^2^^^ + 0(7.-1), ^ ^ 



where 



.o-^arctan^ 



4 Conclusion 

In this paper we have extended results in [11] in several directions. First of all, for a very large class 
of Levy processes (having infinitely many parameters) we have proved that the Wiener-Hopf factors 
can be expressed as infinite products of linear factors and that the distribution of Se(q) and Ie(q) can 
be identified as an infinite mixture of exponential distributions. Second, we have introduced five eight- 
parameter families of Levy processes, which have a wide range of behavior of small jumps, including one 
family having jumps of finite activitity (the density of the Levy measure has a singularity at zero of order 
X = |), two families with jumps of infinite activity but finite variation (x = 1 or x = |) and two families 
with jumps of infinite variation (x = 2 or x = |). We have also derived precise asymptotic expressions 
for the large solutions to (f){z) = q, which are very useful for numerical computations. These five families 
of processes have rather simple form of the characteristic exponent, especially when x is a half-integer, in 
which case "^{z) is given in terms of elementary trigonometric functions. This fact and the availability of 
efficient numerical schemes for computing the Wiener-Hopf factors and the distribution of extrema make 
these processes very interesting for Mathematical Modeling, in particular in the areas of Mathematical 
Finance and Insurance Mathematics. 
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